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ABSTRACT 

Using photometric and redshift data for the Virgo and Coma clusters, we present evidence 
for a correlation between the light concentration of elliptical galaxies (including dwarf ellipticals) 
and the local 3-D (i.e. non-projected) density of the clusters: more concentrated ellipticals are 
located in denser regions. The null hypothesis (i.e. the absence of any relation) is rejected at 
a significance level better than 99.9%. In order to explain the observed relation, a power law 
relating the galaxy light concentration and the cluster 3-D density is proposed. We study how 
the projection effects affect the form and dispersion of the data-points in the light concentration- 
projected density diagram. The agreement between our model and the observed data suggests 
that there is a paucity of dwarf elliptical galaxies in the cluster central regions. 
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1. Introduction 

Analyzing the Abell Cluster 2443, Trujillo et al. 
(2001, hereafter T01) showed evidence for a cor- 
relation between galaxy light concentration and 
local cluster surface density for elliptical galaxies: 
more centrally concentrated ellipticals appear to 
populate denser regions. If this relation is shown 
to hold for other clusters, it means that the qual- 
itative morphology density relation noted by pre- 
vious authors (e.g. Dressier 1980; Dressier et al. 
1997; Fasano et al. 2000) can be placed on a more 
quantitative basis, that is, the detailed structure 
of the individual galaxies (beyond the broad ellip- 
tical/spiral distinction) is related to their immedi- 
ate environment /density. 

That the structural properties of ellipticals are 
related to the properties of their parent clusters 
was noted early on by Strom & Strom (1978). 
The characteristic sizes of galaxies decrease by a 
factor of 1.5 in the denser regions. This effect 
was explained by tidal disruption and high-speed 
impulse encounters. However, as we discussed in 



T01, this mechanism does not seem to be the cor- 
rect explanation for the correlation presented in 
our previous paper. Since mergers tend to increase 
the concentration of galaxies (White 1983; Barnes 
1990, 1992) and bulges (Aguerri et al. 2001), we 
proposed this mechanism as a possibly more viable 
explanation. 

By understanding how the properties of galax- 
ies relate to those of their parent clusters, we can 
hope to learn about the formation and evolution 
of both. In this paper, the connection between 
galaxy light concentration and cluster density is 
examined and confirmed for the Virgo and Coma 
clusters. We pay special attention to the effects 
of projection on this relation, and show how its 
shape and its scatter can be easily explained. 

2. Galaxy data and measurements 

The data for the Coma cluster were taken from 
a quantitative morphological analysis of the galax- 
ies placed in the central region of this cluster, cov- 
ering an area of 0.28 square degrees (see details 
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in Gutierrez et al. 2002). The images were ob- 
tained on 2000 April 25 and 27 using the Wide 
Field Camera at the 2.5 m Isaac Newton Telescope 
at the Observatorio del Roque de los Muchachos 
on La Palma. The pixel scale is 0".333 pixel -1 , 
and the seeing was I'M. The field was observed 
through R filter for a total integration time of 3900 
s. B-R color information for the observed galaxies 
was obtained from the Coma cluster galaxies cata- 
log presented in Godwin, Metcalfe & Peach (1983). 
This catalog is complete down to tub — 20. 

Briefly, the structural parameters of the ellipti- 
cal galaxies were obtaining by fitting a 2D Sersic 
model r 1 /™ (Sersic 1968) to the observed galax- 
ies. Both the ellipticity shapes of the galaxies and 
the effects of seeing on the images were taken into 
account when fitting the model (details of the pa- 
rameter recovering method are explained in T01). 
We used a Moffat function with (3=2.5 to describe 
the point-spread function. The parameters were 
estimated with an error less than 10% down to 
R=17 (itib = 19), which we assume as our limiting 
magnitude for an accurate morphological struc- 
ture analysis. Assuming H = 75 km s -1 Mpc -1 
(which we do throughout), and a redshift for Coma 
z = 0.023, this implies an absolute limiting magni- 
tude Mb = —15.84 and an observed field covering 
the inner 500 kpc. 

We then computed the concentration index of 
the best-fitting r 1 /™ models using a new index pre- 
sented in Trujillo, Graham & Caon (2001) and fur- 
ther developed in Graham, Trujillo & Caon (2001). 
This index measures the light concentration within 
a profile's half-light radius (r e ): it is the ratio of 
flux inside some fraction a of the half-light radius 
to the total flux inside the half-light radius. For an 
r V™ model, this index can be analytically defined 
as 

7 (2n,6„a 1 /») 



Cr» = 



(1) 



7(2n,&„) 

where n is the shape parameter of the r 1 /™ model 
and b n is derived numerically from the expres- 
sion T(2n) = 27(2n, b n ), with T(a) and -f(a, x) re- 
spectively the gamma function and the incomplete 
gamma function (Abramowitz & Stegun 1964). 
The parameter a can be any value between and 
1, and defines what level of concentration is being 
measured. We used a value of a = 0.3. 

We computed the local cluster surface density 



around each galaxy in our sample in the follow- 
ing way. Using the position information from the 
catalog by Godwin ct al. (1983), we compute the 
distance to the 10th nearest neighbor, rio, and de- 
rived the density as p pro j = 10/ (it x rf ). In order 
to avoid contamination from field galaxies and to 
assure uniform completeness on the whole area, 
we selected from that catalog only those galax- 
ies with m B < 20 (M B < -14.84), and satisfy- 
ing the redshift condition 4000 < cz < 10000 km 
s -1 (a 3a interval around the mean cluster red- 
shift). The redshift information was obtained from 
M. Colles (private communication) from the data 
used in Edwards et al. (2002). For those galaxies 
lacking velocity data, we use the color constraint 
1 < B - R < 2 (see e.g. Mobasher et al. 2001). 

The same calculations were carried out for the 
Virgo Cluster. The concentration indexes were 
computed from the best-fitting Sersic n's pub- 
lished in Caon, Capaccioli & D'Onofrio (1993), 
and in Binggeli & Jerjen (1998); only galaxies 
classified as ellipticals or dwarf ellipticals brighter 
than ttib = 15.45 were used. The projected den- 
sities were computed using the galaxies listed in 
the Virgo Cluster Catalog (Binggeli, Sandage & 
Tammann 1985), from which we selected only 
those with confirmed membership and brighter 
than ms = 16.45. The cut-off magnitudes used 
in Virgo correspond to the same limiting ab- 
solute magnitudes used for Coma (assuming a 
Coma- Virgo distance modulus of 3.50 mag from 
D'Onofrio et al. 1997). 

3. A CV c (l/3) p(3D) power law relation 

The relation between galaxy light concentra- 
tion and local cluster surface density is shown 
on Fig. 1. Solid points represent those ellipti- 
cal galaxies with Mb < —17.5 (22 from Virgo 
and 15 from Coma) whereas dwarf galaxies are 
denoted by open points (21 from Virgo and 42 
from Coma). Since CV e (l/3), as defined in equa- 
tion 1, is a monotonic function of the global shape 
parameter n (Trujillo et al. 2001a), the local den- 
sity - C re (l/3) relation implies a relation between 
the local density and n as well. The strength of 
these correlations are of course equal. We prefer 
to maintain the discussion in terms of the concen- 
tration parameter because it has a more tangible 
meaning that the index n. 
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A point we want to emphasize is the "trian- 
gular" form of the shape observed between the 
galaxy light concentration and the cluster pro- 
jected density. This kind of relation is what one 
would expect in case a relation between galaxy 
light concentration and cluster 3-D density is as- 
sumed. To understand this, we note that projec- 
tion effects will tend to mix both more and less 
concentrated galaxies at higher projected densi- 
ties, but at lower projected densities only galaxies 
with low concentrations will be seen, i.e. no high- 
concentration objects will appear in low density 
environments. 

To illustrate this, we study what the shape and 
the scatter in the correlation look like by con- 
structing 3-D simulated spherical clusters. For 
simplicity, we use the generalized King model den- 
sity profile for the 3-D galaxy distribution: 

P{r) = 2 P+1/2 (2) 

p(r c ) (1 + (r/r c ) 2 )^ +1/2 > 

where r c is the core radius and (3 is a parameter 
that models the "tail" of the profile. The extension 
of the tail decreases with increasing (3. This model 
has a flat (core) behavior at r < r c . 

Each galaxy in the realization has associated 
a galaxy light concentration as a function of the 
cluster 3-D density as follows: 

CrMr)/p(r c )) = a,™^^) (3) 

where C rctTnax and M are the maximum values of 
the concentration and the cluster density, respec- 
tively The maximum value of C re is reached at 
the center r = 0. Based on the highest values of 
the concentration that we observe in our measure- 
ments, we set Cr c ,max = 0.7. On the other hand, 
M=2 /3+i/2 for thc King mo del. 

In our models we assume that the distribution 
of the artificial galaxies which are used to evalu- 
ate the local density (i.e. thc deeper sample) is 
described also by thc King model. To calibrate 
the density of this deep sample we have imposed 
to the model that logp pro j(0)=3 (where p pro j(0) is 
the projected density at the center of the cluster). 
This value is close to the highest density that we 
observe in our observational data (see Fig. 1). 



Following Eq. (3), each point (log(p pro j), C r J 
in Fig. 1 can be understood (one-to-one) in terms 
of the pair (R,r), where r is the 3-D radial dis- 
tance and R is the projected 2-D radial distance 
associated to this point in the cluster model. 

The comparison between the simulated model 
point distribution and the observed point dis- 
tribution is done by using thc generalization 
of the Kolmogorov-Smirnov (K-S) test to two- 
dimensional distributions (Fasano & Franceschini 
1987). Thc two-dimensional K-S statistic used 
(D) is defined as the maximum difference (rang- 
ing both over data points and over quadrants) of 
the corresponding fractions between the data and 
the model. In order to match the observed point 
distribution intervals (see Fig. 1), we restrict the 
above comparison to the 0.18 < C Te < 0.70 and 
1.4 < log/9 pro j < 3 intervals. 

The cluster simulations - described in the Ap- 
pendix - are repeated 1000 times, so as to build 
a distribution function DF for the quantity Z? max 
(the largest difference of D for each cluster simula- 
tion). Finally, we repeat the process, this time us- 
ing our observed data-points, to obtain D maXi0 b s - 

Following standard statistical methods, we 
evaluate the probability associated with the mea- 
sured value of D maX;0 bs by calculating the frac- 
tional area under the DF curve where the fre- 
quency is equal to or less than that of -D m ax,obs- 

On Fig. 2 we plot the isocontours of confidence 
level 68% and 95% associated to our models. To 
illustrate how well one of the "acceptable" models 
can reproduce thc shape of the data distribution, a 
realization of the model (3 — 1.3 and 5 — 0.2 (using 
1000 points) is over-plotted on the observed dis- 
tribution in the logp pro j — C re diagram in Fig 1. 
It clearly shows that the shape and the dispersion 
present in Fig. 1 is just the product of the pro- 
jection effects, even when starting from perfectly 
noiseless 3-D relations, such as the one in Eq. 3. 

By looking at Fig. 1, the toy galaxies crowd 
near the diagonal (in particular the upper right 
corner), while most of the real galaxies (in particu- 
lar dEs) are located in the bottom part of the plot. 
To explain this different distribution we must note 
first that when constructing thc artificial clusters 
we ignore the actual luminosity function of Es and 
dEs in the cluster (i.e. artificial galaxies are placed 
in the clusters just taking into account the density 
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Fig. 1. — The light concentration index is plotted 
versus the projected surface density. Solid points 
correspond to bright ellipticals (Mb < —17.5) in 
the Coma Cluster (triangles) and Virgo Cluster 
(circles) . Open points are dwarf galaxies in Coma 
(triangles) and Virgo (circles). Overplotted on the 
observed points is a simulated cluster realization 
of 1000 points (small dots) following Eq. 3, with 
(3 = 1.3 and 8 = 0.2 (see text for details). 
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Fig. 2. The contours of confidence level 68% 
(solid line) and 95% (dashed line) associated to 
the proposed model are shown. 



profile). 

It is clear that any improved version of the 
present models (out of the scope of this paper) 
would have to take into account the fact that the 
dE galaxies are more abundant than the elliptical 
galaxies. Moreover, as the dE galaxies are also less 
concentrated objects than the E galaxies (see e.g. 
the luminosity-concentration diagram in Fig. 8 of 
Graham et al. 2001), one can expect that in ob- 
served galaxy distributions as the one presented 
in Fig. 1 most of the galaxies are found in the 
bottom part of the plot. 

The distribution of artificial galaxies is also a 
function of the assumed symmetry used to de- 
scribe the cluster (in our case spherical symme- 
try). The diagonal line in Fig. 1 is populated by 
the artificial galaxies which lie over the perpen- 
dicular plane to the line of sight passing through 
the cluster center. The fact that the structure of 
real clusters departs significantly from the spher- 
ical simmctry of our simple model accounts for 
the differences between the observed and simu- 
lated distributions. 

To evaluate the null hypothesis (i.e. the absence 
of any relation between the galaxy light concentra- 
tion and the cluster density) we have constructed 
3-D clusters where each galaxy has a concentra- 
tion index independent of the density and ran- 
domly distributed between < C Te < 0.7. The 
null hypothesis is rejected at a confidence level 
>99.99%. This is the main result of this paper: 
whatever the exact form of the relation between 
galaxy light concentration and cluster density may 
be, it is clear that such a relation exists. Also, 
we have checked the consistency of our data by 
making an internal comparison (based on a 2-D 
K-S test) between the distribution of Virgo and 
Coma galaxies. The hypothesis that the two sam- 
ples come from the same parent population turns 
out to be acceptable (rejection probability ~0.3). 

4. Discussion 

A discussion of the mechanisms that may be 
acting to generate the above correlation was pre- 
sented in T01. There, we explained that a mech- 
anism based on mergers is favored over other sce- 
narios based on tidal friction and high-speed im- 
pulse encounters. In any case, in the center of 
galaxy clusters different mechanisms may be op- 
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erating to reshape the form of the galaxies, and 
no single mechanism is expected to account for 
the whole range of observed properties. 

From the model that we have proposed to anal- 
yse the distribution of galaxies in Fig. 1, it follows 
that dE galaxies do not populate the central region 
of the clusters. Using Eqs. 2 and 3 it is easy to 
determine the radial distance to the center of the 
cluster once the galaxy concentration is known. 
Assuming that the maximum galaxy light concen- 
tration of dEs is 0.4 (see. Fig. 1) and that the 
acceptable models are in the range 1.2 < /3 < 1.8 
and 0.15 < S < 0.35 (see Fig. 2), dE galaxies are 
removed from the center of the clusters out to a 
radius <~ 1 — 3r c . 

Other authors (Seeker, Harris & Plummer 1997; 
Gregg & West 1998; Adami et al. 2001; Andrcon 
2002) have found independent evidence which sup- 
ports this result. This would be an indication that 
in the denser cluster environments the dEs (which 
are only weakly gravitationally bound) can be de- 
stroyed by tidal disruption. For the King model, 
we can make a direct comparison between our re- 
sults and the cluster evolutionary model proposed 
by Merritt (1984), in which tidal forces in the clus- 
ter center disrupt the galaxies. The minimum 
size that a galaxy must have in order not to be 
tidally disrupted is ~ 15/j^ 1 kpc. The maximum 
of the tidal disruption forces is reached at r <~ r c . 
Thus, this scenario predicts that dEs are easily 
destroyed, while the larger and more massive el- 
liptical galaxies are able to survive in the center 
of the cluster. Our results are in good agreement 
with this. 

Another important result which follows from 
the kind of law suggested in Eq. 3 is the possibil- 
ity of making a 3-D reconstruction of the elliptical 
galaxies in the cluster. If the spherical assump- 
tion for the density profile of the cluster is good 
enough, then by measuring the galaxy light con- 
centration of each elliptical it will be possible to 
determine its 3-D radial distance. Of course, this 
kind of reconstruction is a priori possible for any 
law which relates a galaxy property (e.g. concen- 
tration, size, color, etc.) with the cluster environ- 
ment density. 
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versions of this manuscript. Finally, the authors 
are grateful to the anonymous referee for the valu- 
able comments that helped us to improve the rigor 
and clarity of this paper. 

REFERENCES 

Abramowitz M., Stegun I., 1964, Handbook of 
Mathematical Functions. Dover, New York 

Adami, C, Mazure, A., Ulmcr, M.P. & Savine, C, 
2001, A&A, 371, 11 

Aguerri, J. A. L., Balcells, M. & Peletier, R. F., 
2001, A&A, 367, 428 

Andreon, S., 2002, A&A, 382, 821 

Barnes, J. 1990, in Dynamics and Interactions of 
Galaxies, ed. R. Wielen (Heidelberg: Springer- 
Verlag), 232 

Barnes, J. 1992, ApJ, 393, 484 

Binggeli, B., Sandagc, A., Tammann, G. A., 1985, 
AJ, 90, 1681. 

Binggeli, B. & Jerjen, H. 1998, A&A, 333, 17. 

Caon, N., Capaccioli, M., & D'Onofrio, M. 1993, 
MNRAS, 265, 1013 

D'Onofrio, M., Capaccioli, M., Zaggia, S.R., & 
Caon, N. 1997, MNRAS, 289, 847 

Dressier, A. 1980, ApJ, 236, 351 

Dressier, A., Ocmlcr, A. Jr., Couch, W. J., Smail, 
I.,Ellis, R. S., Barger, A., Butcher, H., Pog- 
gianti, B. M., Sharpies, R. M. 1997 ApJ, 490, 
577 

Edwards, S. A., Colless, M., Bridges, T J., Carter, 
D., Mobashcr, B., Poggianti, B. M., 2002, 
ApJ,567, 178. 

Fasano, G., & Franceschini, A., 1987, MNRAS, 
225, 155 



5 



Fasano, G., Poggianti, B. M., Couch, W.J., Bet- 
toni, D., KJaergaard, P., & Moles, M. 2000, 
ApJ, 542, 673 

Godwin, J. G., Metcalfe, N., Peach, J. V., 1983, 
MNRAS, 202, 113. 

Graham, A. W., Trujillo, L, & Caon, N., 2001, AJ, 
122, 1707 

Gregg, M. D., West, M. J. 1998, Nature, 396, 549 

Gutierrez, CM., Trujillo, I., Aguerri, J.A.L., Col- 
less, M., Caon, N., Cepa, J., 2002, in prepara- 
tion 

Merritt, D., 1984, ApJ, 276, 26 

Mobasher, B. et al, 2001, ApJS, 137, 279 

Seeker, J., Harris, W. E., Plummcr, J. D., 1997, 
PASP, 109, 1377 

Sersic, J. 1968, Atlas de Galaxias Australcs 
Cordoba: Obs. Astronomico 

Strom, S. E. & Strom, KM. 1978, ApJ, 225, L93 

Trujillo, I., Graham A.W. & Caon, N., 2001a, MN- 
RAS, 326, 869. 

Trujillo, I., Aguerri, J. A. L., Gutierrez, C. M. & 
Cepa, J. 2001, AJ, 122, 38 (T01) 

White, S.D.M. 1983, in Dynamics and Interactions 
of Galaxies, ed. E. Athanassoula (Dordrecht; 
Reidel), 337 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 



6 



A. The cluster simulations 

We carried out cluster simulations as follows. 

1. We populate each cluster realization with the same number of artificial galaxies as the observed ones 
(99), having radial distances in accord with the density law in Eq. 2, and assign to them the concen- 
tration index given by Eq. 3. The data-points are then projected on the log /9 pr oj — C re diagram. 

2. For each point (xi,yi) in the diagram 1 , we compute the fractional number of data-points in the four 
quadrants f a (x > x u y > y t ), f b (x <x u y> y t ), f c (x < x u y < yi) and f d (x > x,,y < y,). 2 . 

3. We then compute, for the same points, the expected fractions given by our density model, which can 
be determined analytically using the expressions: 

f _ jp, «x Sr' cos(&icsm(R* /u))p(u)u 2 du 
Jo P(u)u 2 du 

f a + f b = F(r*) (A2) 

fc + fd = l-F(r*) (A3) 

, f°? cos(a,rcsm(R* /u))p(u)u 2 du 

f d = l-F(r*)-± ) -LJp (A4) 

J p(u)u 2 du 

where (R*, r*) are the values of projected and the 3-D radial distances respectively associated to the 
observed point (log p pTO j (i) , C re (i)), and F(r) is the cumulative distribution function defined (for a 
radial symmetric model) as: 

F(r) = % \ \ 2 . A5 
J p{u)u 2 du 

with p(u) given by Eq 2. For the King model proposed above, F(r) can be derived analytically when 

l</3<2: 

p( r \ _ 4 /rV ^q^ + ^F^sA/j + iAsA-^) 2 ) 
t[r} 3V^\r c ) r(/3-i) 

where 2F1 {a, b; c, z) is the hypergeometric function (Abramowitz & Stegun 1964, p. 556). 

4. For each point (x i} yi) we find the maximum difference between the observed and expected fractions 
in the four quadrants. This is repeated for all data-points, so as to obtain the largest difference -D max 
for each cluster simulation. 



1 In our case (log p pro j («), C re (i)). 

2 By construction, f a + fb + fc + fd = 1- 
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